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Abstract 

Bertrand's theorem asserts that any spherically symmetric natural Hamiltonian sys- 
tem in Euclidean 3-space which possesses stable circular orbits and whose bounded tra- 
jectories are all periodic is either a harmonic oscillator or a Kepler system. In this paper 
we extend this classical result to curved spaces by proving that any Hamiltonian on a 
. spherically symmetric Riemannian 3-manifold which satisfies the same conditions as in 

Bertrand's theorem is superintegrable and given by an intrinsic oscillator or Kepler sys- 
tern. As a byproduct we obtain a wide panoply of new superintegrable Hamiltonian 
systems. The demonstration relies on Perlick's classification of Bertrand spacetimes and 
on the construction of a suitable, globally defined generalization of the Runge-Lenz vec- 
tor. 

PACS: 02.30.1k, 02.40.Ky, 45.20.Jj 

Keywords: Integrable Hamiltonian systems, Runge-Lenz vector, spherical symmetry, 
static spacetimes, periodic orbits, Kepler problem. 

q ; 1 Introduction and preliminary definitions 

The Kepler problem and the harmonic oscillator are probably the most thoroughly studied 
systems in classical mechanics. The reasons for this are twofold. First, these potentials play 
a preponderant role in Physics due their connection with planetary motion and oscillations 
^ . around a nondegenerate equilibrium. Second, these potentials are of particular mathemati- 

cal interest due to the existence of additional (or "hidden") symmetries yielding additional 
constants of motion. In fact, both the Kepler and the harmonic oscillator Hamiltonians 
are (maximally) superintegrable in the sense that they have the maximum number (four) of 
functionally independent first integrals other than the Hamiltonian0 

Bertrand's theorem [6] is a landmark result which characterizes the Kepler and harmonic 
oscillator Hamiltonians in terms of their qualitative dynamics. A precise statement of this 
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theorem is given below. We recall [18] that the first condition, which is occasionally forgotten, 
is necessary in order to exclude potentials of the form V(q) = — If ||q||~ s , with K > and 
s = 2,3,... 

Theorem 1 (Bertrand). Let H = ^||p|| 2 + V(q) be a natural, spherically symmetric Hamil- 
tonian system in a domain o/IR 3 . Let us suppose that: 

(i) There exist stable circular orbits, 
(ii) All the bounded trajectories are closed. 

Then the potential is either a Kepler (V(q) = -A/||q|| + B) or a harmonic oscillator potential 
(V(q) = ^4||q|| 2 +B). In particular, H is superintegrable. 

Analogues of the Kepler and harmonic oscillator systems in curved spaces have been of 
interest since the discovery of non-Euclidean geometry. In fact [49] , the "intrinsic" Kepler and 
harmonic oscillator problems on spaces of constant curvature were studied by Lipschitz and 
Killing already in the 19th century, and later rediscovered by Schrodinger [48] and Higgs [25]. 
In both cases it was established that these systems are superintegrable and satisfy Properties 
(i) and (ii) above. 

A considerably more ambitious development was Perlick's introduction and classification 
of Bertrand spacetimes [47], which was based on the following observation. Let (M,g) be 
a Riemannian 3-manifold and consider the space A4 = M x R endowed with the warped 
Lorentzian metric rj = g — ydt 2 , with V a smooth positive function on M. Then the trajec- 
tories in (A4,rj), that is, the projections of inextendible timelike geodesies to a constant time 
leaf M x {to}) correspond to integral curves of the Hamiltonian H = ^\\p\\^ + V(q) in (the 
cotangent bundle of) M. Thus Perlick introduced the following 

Definition 2. A Lorentzian 4-manifold (M x R, rj) is a Bertrand spacetime if: 

(i) It is spherically symmetric and static in the sense that r) = g — ydt 2 and M is dif- 
feomorphic to {r\,T2) x 8 2 , where the smooth function V depends only on r and the 
Riemannian metric g on M takes the form 

g = h(r) 2 dr 2 + r 2 (d9 2 + sin 2 9 d(p 2 ) (1) 

in the adapted coordinate system (r, 6, ip). Here r\,T2 € K + U {+00} . 

(ii) There is a circular (r = const.) trajectory passing through each point of M. 

(iii) The above circular trajectories are stable, that is, any initial condition sufficiently close 
to that of a circular trajectory gives a periodic trajectory. 

Perlick's main result was the classification of all Bertrand spacetimes, recovering the 
classical Bertrand theorem as a subcase. However, two main related questions remained to 
be settled. On the one hand, the potentials V in Perlick's classification lacked any physical 
interpretation, and this was in strong contrast with the Euclidean case. This drawback 
was circumvented in Ref. [3], where we showed that the two families of Perlick's potentials 
correspond to either the "intrinsic" Kepler or harmonic oscillator potentials in the underlying 
3-manifold (M,g). On the other hand, the issue of whether the corresponding Hamiltonian 
systems were superintegrable in some reasonable sense was left wide open. In fact, Perlick's 
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only remark in this direction was that, by virtue of a theorem of Hauser and Malhiot [24], 
only two concrete models among the family of Bertrand spacetimes admitted a quadratic 
additional integral coming from a second rank Killing tensor. 

A careful analysis of the literature reveals that many particular cases of Bertrand metrics 
have been thoroughly analyzed and shown to be superintegrable [21,22,27,28], and that in 
many cases they have been shown to admit a generalization of the classical Runge-Lenz vec- 
tor as an additional first integral. The physical and mathematical interest of these models 
(and thus of Bertrand spacetimes) is fostered by their connections with the theory of mag- 
netic monopoles, with differential and algebraic geometry, and with low-dimensional mani- 
fold theory [9,10,33-35,43,45,46,50]. The relation between Bertrand spaces and monopole 
motion is not totally incidental. Indeed, an ample subclass of Bertrand spacetimes admit- 
ting some kind of generalized Runge-Lenz vectors (the so-called multifold Kepler systems) 
were introduced by Iwai and Katayama [27, 28] as generalizations of the Taub-NUT met- 
ric, whose geodesies asymptotically describe the relative motion of two monopoles (see, for 
instance, [1,7,11,29,39,40]). Interestingly, superintegrable Hamiltonian systems on curved 
spaces have recently attracted considerable attention also within the integrable systems com- 
munity, especially in low dimensions (cf. [2,4,5,30-32] and references therein). 

The main result of this article is that all Bertrand spacetimes are indeed superintegrable, 
their superintegrability being linked to the existence of a generalized Runge-Lenz vector. 
This enables us to present an optimal version of Bertrand's theorem (Theorem 1 16j) on spheri- 
cally symmetric manifolds which includes the classification of the natural Hamiltonians whose 
bounded orbits are all periodic [47], the physical interpretation of the corresponding potentials 
as Kepler or harmonic oscillator potentials, in each case, and the proof of the superintegra- 
bility of these models. This settles in a quite satisfactory way a problem with a large body 
of previous partial results scattered in the literature. 

It is standard that the superintegrability of the Kepler system stems from the existence 
a conserved Runge-Lenz vector, whose geometric significance is described from a modern 
perspective in [23]. On the other hand, the superintegrability of the harmonic oscillator is 
usually established either using explicit (scalar) first integrals or the conserved rank 2 tensor 
C = 2cj 2 q (8) q + p ® p, which is sometimes preferable for algebraic reasons [19]. That the 
latter approach is closely related to a (multivalued) analogue of the Runge-Lenz vector was 
firmly established in [26]. Motivated by this connection, we have based our approach to 
the integrability of the Bertrand systems on the construction of a generalized Runge-Lenz 
vector, globally defined on a finite cover of M. This construction relies on a detailed analysis 
of the integral curves of the appropriate Hamiltonians. The literature on generalizations of 
the Runge-Lenz vector for central potentials on Euclidean space is vast (see the survey [36] 
and references therein), but unfortunately several interesting papers are severely flawed by 
the lack of distinction between local, semi-global and global existence. 

The article is organized as follows. In Section [2] we recall the two families of Bertrand 
spacetimes entering Perlick's classification, which are labeled by two coprime positive inte- 
gers n and m. We also include the characterization of Perlick's potentials as the intrinsic 
Kepler or harmonic oscillator potentials of the corresponding Riemannian 3-manifolds (M, g) 
and briefly discuss several physically relevant examples. In Section [3] we consider the asso- 
ciated natural Hamiltonian systems on (M,g) and compute their integral curves in closed 
form (Proposition [7|). Using this result we easily derive that the latter Hamiltonians are geo- 
metrically superintegrable (cf. Definition [9] and Proposition fTUT) in the region of phase space 
foliated by bounded orbits, as happens with the harmonic oscillator and Kepler potentials in 
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M 3 . Our central result is a stronger superintegrability theorem (Theorem I12p that we present 
in Section 01 where we construct a generalized Runge-Lenz vector globally defined on an 
n-fold cover of M. As a corollary of this construction we also obtain a global rank n tensor 
field in M invariant under the flow and a wide panoply of new superintegrable Hamiltonian 
systems. Lastly, in Section [5] we combine the results established in the previous sections to 
obtain an optimal extension of Bertrand's theorem to curved spaces (Theorem I16p, 



2 Harmonic oscillators and Kepler potentials in Bertrand space- 
times 

In this section we shall define the "intrinsic" Kepler and harmonic oscillator potentials in 
a spherically symmetric 3-manifold and show how Bertrand spacetimes are related to the 
Kepler and harmonic oscillator potentials of any of its constant time leaves. Most of the 
material included here is essentially taken from Ref. [3]; for the sake of completeness, let us 
mention that further information on geometric properties of Green functions can be consulted 
e.g. in [13-15,37,38]. 

We start by letting (M, g) be a Riemannian 3-manifold as in Definition [2J In particular, 
the metric g takes the form 

ds 2 = h(r) 2 dr 2 + r 2 (d6 2 + sin 2 6 d^ 2 ) . (2) 

It is standard that if u(r) is function which depends only on the radial coordinate, then its 
Laplacian is also radial and reads as 

1 d / r 2 du 

A g u(r) 



r 2 h(r) dr \ h(r) dr J 

As the Kepler potential in Euclidean three-dimensional space is simply the radial Green 
function of the Laplacian and the harmonic oscillator is its inverse square, it is natural to 
make the following 

Definition 3. The (intrinsic) Kepler and the harmonic oscillator potentials in (M,g) are 
respectively given by the radial functions 

V K (r) = A l(J a r'- 2 h(r') dr' + B^j , V R (r) = A 2 (J r'- 2 h(r') dr' + B^j , (3) 
where a,Aj,Bj are constants. 

Example 4. Let (M,g) be the simple connected, three-dimensional space form of sectional 
curvature n. In this case the metric has the form ^ with 

h{r)° ' 



1 — Kr 2 

The corresponding Kepler and harmonic oscillator potentials are therefore 

V K = Vr~ 2 -K, Vu = —^ (4) 



r 2 — k 
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up to additive and multiplicative constants. In terms of the distance function p K to the point 
r = this can be rewritten as 

t/ r + ( r \ \r tan 2 (^p K ) 

V K = COt [y/K p K ) , K H = , 

thus reproducing the known prescriptions for the sphere and the hyperbolic space [5,49]. The 
Euclidean case is recovered by letting k — > 0. 

Now let us consider the spherically symmetric spaces (M,gj) (j = I, II) defined by the 
metrics 

Type I: ds 2 = — + r 2 {d9 2 + sin 2 dc^ 2 ) (5a) 

n z (1 + Ar) 



2m 2 1 - Dr 2 ± ^(1 - Dr 2 ) 2 - Kr 4 
Type II: d, 2 = V n 2 ((1 - Dr 2 ) 2 - Kr*) ^ + r 2 (d0 2 + sin 2 #d^) (5b) 

where D and iT are real constants and m and n are coprime positive integers. The maximal 
interval (ri,r2) can be easily found from these expressions. These Riemannian 3-manifolds, 
which first appeared in [47] (where the quotient n/m was called (3), will be henceforth called 
Bertrand spaces. A short computation shows that, up to a multiplicative constant, the Kepler 
potential of a Bertrand space of type I is 



Vi = \fr- 2 + K + G , (6a) 

whereas the harmonic oscillator potential of one of type II can be written in the convenient 
form 

V n = G =p r 2 (l - Dr 2 ± - Dr 2 ) 2 - KrA . (6b) 

Here G is an arbitrary constant. 

By comparing with Ref. [47], the above digression immediately yields the following 

Proposition 5. (Ai,rj) is a Bertrand spacetime if and only if it is isometric to the warped 
product (M x M,gj — %r), with (M,gj) a Bertrand space of type j (j = I, II, cf. ©J and Vj 
given by ([6]). 

In particular, this shows that Per lick's obtention of two different kinds of Bertrand space- 
times has a natural interpretation [3]: they are associated to either Kepler (type I) or har- 
monic oscillator (type II) potentials. The multiplicative constant of the potentials is inessen- 
tial and can be eliminated by rescaling the time variable. 

Example 6. We conclude this section with a brief discussion of a few examples of physically 
relevant spaces that are Bertrand. This intends both to serve as motivation and to help the 
reader gain some insight on Bertrand spaces. A more detailed discussion can be found in [3]. 

(i) Spaces of constant curvature. The metric of the simply connected Riemannian 3- 
manifold of constant sectional curvature n is usually written as 

ds 2 = dr 9 + r 2 (d# 2 + sin 2 9 dp 2 ) . 
1 — nr z 
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We have already seen that the Kepler and harmonic oscillator potentials in these spaces 
are given by (j3|), and it is well known that all the bounded integral curves of both 
systems are periodic. This result is immediately recovered by noticing that the Kepler 
system is recovered from the type I Bertrand spacetimes when n = m = 1 and K = —k, 
whereas the harmonic oscillator is obtained as the type II Bertrand spacetime with 
n/m = 2, K = and D = k. 

(ii) Darboux space of type III. Consider the metric 



ds 2 = k2 + 2 ^,t Wk ' ± ±l d r 2 + r 2 (d9 2 + sin 2 6 d^ 2 ) , 
2{k 2 +4r 2 ) v ' 

whose intrinsic harmonic oscillator potential is given by 

2fc 2 r 2 



k 2 + 2r 2 + k\Jk 2 + 4r 2 



up to multiplicative and additive constants. This defines a Bertrand spacetime of type 
II with parameters n/m = 2, K = A/k 4 and D = —2/k 2 . 

Let us introduce coordinates Q = (Q 1 , Q 2 , Q 3 ) as 

/(k 2 + 4r 2 ) 1/2 - £A 1/2 
Q = ( J ( cos 9 cos ip, cos 6 sin ip, sin 0) . 

In terms of these coordinates, the above metric and potential read as 

A; 2 ||Q|| 2 



ds 2 = (k + ||Q|| 2 )||dQ|| 2 , Vn 



fc+||Q| 



Thus we recover the three-dimensional Darboux system of type III [32]. The Dar- 
boux system of type III is the only quadratically superintegrable natural Hamiltonian 
system in a surface of nonconstant curvature which is known to admit quadratically 
superintegrable ^-dimensional generalizations [2]. 

(iii) Multifold Kepler systems. The family of multifold Kepler systems was introduced by 
Iwai and Katayama [27, 28] as Hamiltonian reductions of the geodesic flow in a gener- 
alized Taub-NUT metric. These systems are given by the metrics and potentials 



ds z = ||Q||™^(a + 6||Q||™)||dQ 

Vn = , Tn U 2 ||Q|r 2 + M 2 c||Q||-~ 2 + M 2 d||Q"~~- 



a + 6[|Q| 

with Q = (Q , Q 2 , Q 3 ), a,b,c,d,u constants and n,m coprime positive integers. The 
substitution 

/ (a 2 + 4br 2 )2 - q \ n _ 
(4=1 — I { cos cos ip, cos sin 99, sin tt) 

shows that the multifold Kepler models are equivalent to the type II Bertrand systems 
with parameters K = 4a _4 6 2 and D = — 2%. It should be noticed that the Darboux 
space of type III is a particular case of the multifold Kepler systems. 
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3 The orbit equation and geometric superintegrability 



Hereafter we shall analyze the properties of the Hamiltonian systems in (M,gj) given by 

1 



I, II, 



(7) 



where the metric gj and the potential Vj are respectively defined by (0) and © . As previously 
discussed, the orbits of these systems correspond to trajectories of the associated Bertrand 
spacetimes. It should be noticed that in the adapted coordinate system, these Hamiltonians 
read as 



Hi 



n 



in 



p% 



r- sin 2 i 



H n = - 

2 



((1 - Dr 2 ) 2 - Kr 4 ) p 2 r 



2m 2 1 - Dr 2 ± - Dr 2 ) 2 - Kr 4 



+ Vr~ 2 + K + G. 



r 2 r 2 sin 2 i 



(8a) 



=F r 2 ( 1 - Dr 2 ± y/(l - Dr 2 ) 2 - Kr 4 ) +G 



(8b) 



where p r is the momentum conjugate to r and pe and p v are defined analogously. 

In this section we shall derive the simplest superintegrability property of the Hamiltonian 
systems (J7J (cf. Proposition I10H . which nonetheless seems to have escaped unnoticed so far. 
The proof of this result relies on the fact that, by definition, the orbits of ([7]) define an 
invariant foliation by (topological) circles in an open subset £1 C T*M of the phase space of 
the system. E.g., in the classical Kepler problem 



V= {(q,p) e 



:#(q,p) <0, qxp/0} 



is the set of points with negative energy and nonzero angular momentum, whereas for the 
harmonic oscillator one can take O = (M 3 x M 3 )\{(0, 0)}, i.e., the whole phase space minus 
the equilibrium. In Proposition [7] below we compute the expression of the orbits in closed 
form, revealing that the above foliation is actually a locally trivial fibration. This allows us 
to resort to the geometric theory of superintegrable Hamiltonian systems [12], yielding the 
first superintegrability result for ([7]). 

Before discussing the precise statement of Proposition [TUl let us compute the orbits of 
the Hamiltonian ([7]). In fact, the closed expression that we shall derive is not only used in 
the proof of Proposition [T0| but it is also a key element of Theorem [T2| where a stronger 
superintegrability result is presented. It is convenient to introduce the rectangular coordinates 
q = (q 1 , q 2 , q 3 ) associated to the spherical coordinates (r,0,ip) as 



I r cos a cos ip, r cos a sin ip, r sin 



e). 



(9) 



The conjugate momenta will be denoted by p = {jp\,P2,Pz)- Clearly the coordinates (q, p) 
are globally defined in T*M. We shall use the notation •, x and || • || for the Euclidean inner 
product, cross product and norm in M 3 and call E = Hj(p, q) and J 2 = ||q x p|| 2 the energy 
and angular momentum of an integral curve (q(t),p(i)) of ([7|). Obviously E and J 2 are 
constants of motion. 
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Proposition 7. Let 7 be an inextendible orbit of the Hamiltonian system (JT)) which is con- 
tained in the invariant plane {9 = ?}. Then 7 is given by 



cos ( — — Lfo] = -==k (10a) 



m 



1 + J 2 Vr~ 2 + K 



V / 1 + 2J 2 (E-G) + KJ 4 



if j = I and by 

/nip \ J 2 r~ 2 (l - Dr 2 ± J (I - Dr 2 ) 2 - Kr 4 ) + DJ 2 + 2G-2E , s 

cos ipo = . (10b) 

\m "7 ^/{2E -2G- DJ 2 ) 2 ± 4J 2 - if J 4 

i/j = II. ifere y>o is a real constant. 

Proof. We begin with the case j = I. The crucial observation is that the orbit equation 



2 r\l + Kr 2 )\dip J 1 r 2 



n 

simplifies dramatically with the change of variables 



u = \J r 2 + K , 

in terms of which the potential and the inverse square term read as 

Vi = u + G , r- 2 = u 2 - K . 
The orbit equation is then given by 

m J dU ^ 2 = 2E - 2G + KJ 2 -2u- J 2 u 2 , 



n dip 

which can be readily integrated to yield 



nip \ 1 + J 2 u 

cos 1 ipn = — . 

1 m ™J yJ\ + 2J 2 (E-G) + KJ* 



for some constant ipo. 

When j = II the treatment is analogous. Now the orbit equation reads 

1 - Dr 2 ± - Dr 2 ) 2 - Kr A ( mJ dr \ 2 _ J 2 



DH ± V(l - -Dr 2 ) 2 - & Vm J dr y _ J 2 
r 4 [(l-L>r 2 ) 2 -£> 4 ] \~d^J ~ 11 ~ 2^ ' 



and it is convenient to introduce the variable 



v = r 



- 2 ' 1 - Dr 2 ± ^(l-Dr 2 ) 2 -^ 4 . 



In terms of this new coordinate the potential is simply Vn = whereas the inverse 

square term is given by 

o v 2 + 2Dv + K 

r = . 

2v 

Hence a straightforward computation shows that the orbit equation is 

V^l^L V = 4 (E - G)v - J 2 (v 2 + 2Dv + K)±A, 
n dip J 
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thereby obtaining 



f ncp \ 

ip = 

\m J 



J 2 {v + D) + 2G-2E 



y / (2E -2G- DJ 2 ) 2 ± 4J 2 - if J 4 



Here tpo is a real constant. 



□ 



Remark 8. Eqs. (jlOl) are well denned also when J = 0. Moreover, it is not difficult to check 
that r can be readily expressed as a function of <p by performing some manipulations in the 
right-hand side of (I10p . 

We shall now specify what is understood by geometric superintegrability. Let Fq be a 
smooth Hamiltonian defined on a 2d dimensional symplectic manifold N admitting s > d — 1 
functionally independent first integrals Fi, . . . , F s other than the Hamiltonian. Let us suppose 
that F = (Fq, F%, . . . , F s ) is a submersion onto its image with compact and connected fibers, 
which by Ehresmann's theorem (cf. e.g. [41]) implies that its level sets define a locally trivial 
fibration T of N. If s > d, not all the latter first integrals can Poisson-commute: the usual 
condition to impose is that there exists a matrix-valued function P : F{N) — > Mat(s + 1) of 
rank s — d + 1 such that 



In particular, when s = d — 1 this yields the usual definition of Liouville integrability. Well 
known generalizations of the Liouville-Arnold theorem [42,44] show that every fiber of F is 
an invariant {2d — s — l)-torus, and that the motion on each of these tori is conjugate to a 
linear flow. Moreover, the fibration T has symplectic local trivializations. 

Geometrically, the existence of the function P means that J- has a polar foliation [12], i.e, 
a foliation J 7 ^ whose tangent spaces are symplectically orthogonal to those of T . Similarly, 
the rank condition in Eq. (|lip is tantamount to demand that the invariant (2d — s — l)-tori 
of the foliation be isotropic. Thus the crucial element in the geometric characterization of 
superintegrability is the bifoliation (J 7 , J 7 - 1 ), which is a type of dual pair as defined in [51]. 
One is thus led to introduce the following definition (cf. [12] and the survey [17], where slightly 
different wording is used): 

Definition 9. A Hamiltonian system on a symplectic 2d-dimensional manifold is geomet- 
rically superintegrable with s > d — 1 semiglobal integrals if the Hamiltonian vector field is 
tangent to a locally trivial fibration by isotropic {2d — s — l)-tori which admits a polar folia- 
tion. If s takes the maximum value 2d — 2 we shall simply say that the system is geometricaly 
superintegrable. 

Of course, generally not all the phase space of a (super)integrable system is fibered by 
invariant isotropic tori: there can be, e.g., singular points and unbounded orbits. But it is 
customary and of interest to restrict one's attention to the region where such fibration is well 
defined. In the case when s = d—1 (Liouville integrability), the invariant tori are Langrangian 
and therefore J- 1 - = explaining why the bifibration {J 7 , ^ r_L ) is less well known than the 
fibrations by Lagrangian tori. (However, an advantage of the bifibration is that, under mild 
technical assumptions, it is uniquely determined (and finer), whereas for integrable systems 
with additional integrals there is some arbitrariness in the choice of invariant Lagrangian 
tori.) It should be noticed that the above structure yields "semiglobal" (i.e., defined in a 
tubular neighborhood of each torus) first integrals associated to the existence of generalized 



{F i ,F j } = P ij oF 



< i, j < s. 



(11) 
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action-angle coordinates; a detailed account can be found in [8,12,17]. The content of the 
following proposition is that the Bertrand systems ([7]) are geometrically superintegrable in 
the region foliated by periodic orbits. 

Proposition 10. Let £1 be the region ofT*M where all the orbits of Hj are periodic. Then 
is geometrically superintegrable. 

Proof. It easily follows from Proposition[7Jthat the orbits of Hj define a locally trivial fibration 
by (topological) circles in f2. The fibers are certainly isotropic, as they are one-dimensional, 
and the flow of Hj on each fiber is conjugate to the linear one because Hj does not possess 
any critical points in S7. Moreover, it stems from Proposition [7] that the function — > R + 
mapping each point in f2 to the length of the (periodic) orbit passing through it is smooth, 
which in turn readily implies that the period function is also smooth and nonvanishing in 
this region. Hence a theorem of Fasso [16] implies that Hj is geometrically superintegrable, 
proving the proposition. □ 



4 The generalized Runge— Lenz vector 

In this section we shall prove a stronger superintegrability result for the Bertrand Hamil- 
tonians ([7]). More precisely, we shall provide a semi-explicit construction of an additional 
vector first integral which we shall call the generalized Runge-Lenz vector. This vector field 
is defined on an n-fold cover M of the original space M, and it is invariant under the flow 
generated by the lift of the Bertrand Hamiltonian to the covering space M. In M, this vector 
field induces a global tensor field of rank n which is preserved under the flow of H. As before, 
n is the positive integer which appears in Eq. (J7J). 

As regards the superintegrability properties of the Hamiltonian systems ([7]), the spherical 
symmetry of these systems readily yields three first integrals other than the Hamiltonian, 
which can be identified with the components of the angular momentum. The idea of looking 
for generalizations of the Runge-Lenz vector in order to find an additional integral of motion 
is not new: an updated and rather complete review of the related literature can be found 
in [36]. Here we shall use our information about the integral curves of (JTJ) and some ideas 
already present in the work of Fradkin [20] and Holas and March [26] . 

Let us start by recalling Fradkin's construction [20] of a local vector first integral for the 
Hamiltonian system 

tfo = ^||p|| 2 + t/(||q||), 

where J7([|q[|) is an arbitrary central potential and (q, p) G R 3 x M 3 . The starting point is 
the following trivial remark. Consider an integral curve q(t) of Hq contained in the plane 
{9 = t^} C K 3 , where (r,6,ip) are the usual spherical coordinates. We can assume without 
loss of generality that we have taken the initial condition </?(0) = and use the notation 
r = ||q||, J = p v = r 2 ip. A simple computation shows that the derivative along this integral 
curve of the unit vector field 

cos (p sin 03 , . 

a= -ql-fqxqxp 12 

r rJ 

is identically zero, as in fact a(i) is the constant vector (1,0,0). Fradkin's observation was 
that if cos if and J -1 sin ip can be expressed in terms of q and p in a domain O C M 3 \{0}, then 
the resulting vector field is a first integral of Hq in f2. When Hq is the Kepler Hamiltonian, 
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the generalized Runge-Lenz vector field is well defined globally and essentially coincides with 
the classical Runge-Lenz vector divided by its norm. When Hq is the harmonic oscillator, the 
generalized Runge-Lenz vector is multivalued (this can be neatly understood by considering 
the turning points of the orbits), but can be used to recover the conserved tensor field 
C = 2w 2 q (g> q + p (g) p associated to the SU(3) symmetry [26]. 

Definition 11. Let H be a Hamiltonian system defined on (the cotangent bundle of) a 3- 
manifold N. We say that H admits a generalized Runge-Lenz vector if there exists a nontrivial 
horizontal vector field A in T*N which is constant along the flow of H. 

Obviously the conserved vector A is nontrivial if it is not constant and cannot be written 
in terms of the energy and the angular momentum integrals, and we recall that a horizontal 
vector (resp. tensor) field in T*N can be simply understood as a vector (resp. tensor) field in 
N which depends on both the positions and the momenta. The main problem with Fradkin's 
approach is that, of course, it is not at all obvious how to obtain sufficient conditions ensuring 
that these local integrals are well defined globally, while local sup erintegr ability is trivial in a 
neighborhood of any regular point of the Hamiltonian flow. However, we shall see below that 
Fradkin's approach works well for the kind of Hamiltonian systems that we are considering 
in this paper, and that one can construct a globally defined generalized Runge-Lenz vector 
(cf. Eq. (|16p below) which is roughly analogous to (| 1 2 j) . 

Theorem 12. Consider a Hamiltonian of the form (J7J) ; with m,n coprime positive integers. 
Then there exists an n-fold cover M of M such that the lift of this Hamiltonian to M admits 
a generalized Runge-Lenz vector. 

Proof. We shall call Hj, j = I, II, the Hamiltonian (J7J). Let 7 be an inextendible orbit of Hj, 
which can be assumed to lie in the invariant plane {9 = ^}. By Proposition [71 and taking 
ipo = in Eq. (fT0|) without loss of generality, 7 is the self-intersecting curve given by 

cos^ = X (r 2 ,J 2 ,E), (13) 
m 

where x is the function 

1 + jVr- 2 + K 



X (r 2 ,J 2 ,E) = { 



yfl + 2.P{E -G)+KJ A 



if j = l. 



J 2 r- 2 (l - Dr 2 ± y/(l - Dr 2 ) 2 - Kr 4 ) + DJ 2 + 2G-2E . f . _ ^ 



^/{2E -2G- DJ 2 ) 2 ± 4J 2 - KJ 4 
Moreover, the chain rule immediately yields 

n<p m d ( n(p\ mr d r 2 7 2 tp\ rv • 2 



sin - =-~( cos^ =- — ^ X (r 2 ,J 2 ,E) = Q(rr,r 2 ,J,E), (14) 
m n dtp \ m J n<p or 

where 

2 

Q(rr,r 2 , J, E) = -2rr r ^—(D 1 x)(r 2 ,J 2 ,E) 

nJ 

and Dix stands for the derivative of the function x with respect to its first argument. It 
should be noted that these expressions are well defined also for J = 0. 
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Using the properties of the Chebyshev polynomials it is trivial to express cosny? and 
sin nip in terms of r,f,J and E as 

cos nip = T m (^ cos =T m (x(r 2 ,J 2 ,E)), 

sinnip = sin — C/ m _i( cos — ) = @(rr,r 2 , J, E) U m -i(x(r 2 , J 2 , E)) . 

m \ m J ' 

Here T m and U m respectively stand for the Chebyshev polynomials of the first and second 
kind and degree m. Setting 

S 1 = [ z G C : \z\ = 1} , 
we find it convenient to define the analytic § 1 -valued map 

£ n (rf,r 2 ,J,E)=T m ( X (r 2 ,J 2 ,E))+iG(rf,r 2 ,J,E)U m - 1 {x(.r 2 ,J 2 ,E)) , 
in terms of which the orbit 7 is characterized as 

j n *> = £ n (rr,r 2 ,J,E). (15) 

It stems from Fradkin's argument that (I12p yields a vector first integral of ([8]) in any 
region where e ltp can be unambiguously expressed in terms of the coordinates (q, p). However, 
Eq. (|15|) does not determine the angle ip univocally modulo 2ir because the map z 1— > z n of 
the unit circle onto itself has degree n, so that Fradkin's construction is, a priori, not global. 
As a matter of fact, it is obvious that Eq. (fT5j) only defines <p modulo 27r/n, thus yielding an 
n-valued additional integral. 

The aforementioned problem is a consequence of the fact that the orbit 7 has self- 
intersections. It is standard that this difficulty can be circumvented by means of an appropri- 
ate covering space of our initial manifold. The construction which we shall next outline is in 
fact analogous to that of the Riemann surface of the function z 1— > z n . We denote by j(t) the 
periodic integral curve of ([7]) defined by the orbit 7 C M and take an n-fold cover n : M — > M 
of M such that the lift j(t) of j(t) to M is a smooth path without self-intersections. Notice 
that 7(t) is actually an integral curve of the lifted Hamiltonian 

^' = ^lln* gj + (^° n )(^ .7=1,11, 

where (q,p) E T*M. M is a fiber bundle over M with typical fiber Z n , and for each k G Z n 
we denote by : M — > M the section of M with fiber value k. Obviously A^ is an injective 
map, and an isometry from an open and dense subset M/% C M onto its image in (M,IL*gj). 
One obviously has that n o A^ = id and 

n~\q)= |J A k (q) 

for all q £ M. 

By construction, in each section Afc(M) there exists a determination of the (complex) n-th 
root which allows to solve e 1</3 in terms of e m</3 univocally along Af-^y). Therefore, for each 
k G 7L n there exist real functions Sk and Ck (namely, determinations of arcsin and arccos) 
such that 

e i<p(t) = c k (cosn<p(t)) + iSk(smnip(t)) 
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whenever the point (r(t),9 = ^,<p(t)) lies in A/%(7). Moreover, an easy computation shows 
that the functions 



C k (r 2 ,J 2 ,E) =C k [T m (x(r 2 ,J 2 ,E))) , 
S k (rr, r 2 , J 2 , E) = J- l S k {®{rr, r 2 , J, E) U rn ^ ( X (r 2 , J 2 , E))) 

are analytic in their domains. 

In order to express C k (r 2 , J 2 ,E) and S k (rr, r 2 , J 2 ,E) in a more convenient way, we con- 
sider the lift of the coordinates q to each space A k (M). With a slight abuse of notation, we 
shall still denote these coordinates by q. An immediate computation shows that H* gj\A k (M) 
reads as 

d S 2 = ||dq|| 2 +[/ i (||q||) 2 -l]^^- 2 



where the function h is defined as in Section [21 namely, 

,2 



h{r) 



m 



n 2 {l + Kr 2 ) 

2m 2 (l - Dr 2 ± ^(1 - Dr 2 ) 2 - Kr 4 ^ 
n 2 ((1 - Dr 2 ) 2 - Kr A ) 



if j = I, 



if i = II . 



By differentiation it stems from this formula that the conjugate momentum p to q is given 
by 



p = q-r L"-vii4ii N2 1,1 1 

yielding q = v(q, p) with 



P = q+N|q||) 2 -l]j^q. 



v(q,p)=P+N|q||)- 2 -l]^q. 

As rr = q • q = q • v(q, p), we now have all the ingredients to invoke Fradkin's argument 
(cf. Eq. (I12p . with which (|16|) should be compared) and derive that each component of the 
horizontal vector field A k in T*A k (M) defined by 

A fc = - C fc (||q|| 2 , Hqxpll^tf^p)) q + 5 fc (q-v(q,p), ||q|| 2 , ||q X p|| 2 , ^(q, p)) q X (q X p) 

(16) 

is a constant of motion in A k (M). By construction, the vector fields (with k S Z n ) define 
an analytic global horizontal vector field A in T*M whose Lie derivative along the flow of Hj 
is zero, thereby obtaining the desired unit Runge-Lenz vector. □ 

Remark 13. The particular form of the orbits (|10p and the fact that M is a finite cover of 
M ensure that all the lifted orbits which are bounded are also periodic, and that the lifted 
orbits do not have any self-intersections. Note that if M is endowed with the pulled back 
metric gj = Tl*gj, Hj is a natural Hamiltonian system and II : (M,gj) — > (M,gj) becomes a 
Riemannian cover. 

Corollary 14. Consider a Hamiltonian of the form ([7]) with n = 1. Then the generalized 
Runge-Lenz vector is well defined in all M . 

Proof. It trivially follows from Theorem 1121 □ 
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Corollary 15. Consider a Hamiltonian Hj of the form (J7]), with m,n coprime positive 
integers. Then there exists a horizontal symmetric tensor field in M of rank n which is 
invariant under the flow of Hj . 

Proof. Let us use the same notation as in the proof of Theoremll2l In particular, we consider 
the integral curve 7(i) and the maps used in the proof of Theorem 1121 For each k 6 Z^, 
let us denote by A^(t) the restriction of the horizontal vector field : T*A/%(M) —* R 3 
to the projection of the integral curve j(t) to T*Afc(M). The only observation we need in 
order to prove Corollary [15] is that, by the expression for 7 found in Proposition [7] and the 
definitions of the covering space M and of the horizontal vector fields A^, it easily follows 
that 

A fe (i + |T 7 ) =A k +e{t) 

for all k G Z n , I £ Z, t E R such that 7(f) G M k+i and j(t + ^T 7 ) £ M k . Here T 7 stands for 
the period of the integral curve j(t) and the sum k + £ is to be considered modulo n. This 
periodicity property readily implies that the symmetric tensor product C of Ai, . . . , A n , with 
components 

C^- i ™(q,p)=Af 1 (q,p)---^)(q,p), 

is a well defined, analytic tensor field in M of rank n. As usual, symmetrization of the 
superscripts delimited by curved brackets is understood. To complete the proof of the corol- 
lary, it suffices to notice that C is trivially invariant under the flow of Hj as each A& is a 
(multivalued) first integral. □ 

Some comments may be in order. First, one should observe the dependance of the 
additional integrals (|16p on the momenta is generally complicated (and in particular not 
quadratic), which explains why they are usually so hard to spot [30]. Second, it should be 
noticed that Corollaries [H and [15] yield the usual Runge-Lenz vector and second rank con- 
served tensor (up to a normalization constant) when the Bertrand Hamiltonian we consider 
is the Kepler or harmonic oscillator system in Euclidean space [26,47]. Note, however, that 
given an arbitrary Bertrand Hamiltonian it is usually hard to compute the conserved tensor 
C or the Runge-Lenz vector A in closed form. In this direction, it should be mentioned that 
an additional integral has been explicitly obtained for some of the Bertrand Hamiltonians 
discussed in Example El (cf. e.g. [2,4,21,27] and references therein). 

5 Bertrand's theorem on curved spaces 

In the previous sections we have thoroughly analyzed the superintegrability properties of the 
spherically symmetric natural Hamiltonian systems whose bounded orbits are all periodic. 
When combined with the discussion of harmonic oscillators and Kepler potentials on Bertrand 
spacetimes presented in Section [2] this gives all the ingredients we need to state a fully 
satisfactory analogue of Bertrand's theorem on spherically symmetric spaces: 

Theorem 16. Let H be the Hamiltonian function associated to a Bertrand spacetime, i.e., an 
autonomous, spherically symmetric natural Hamiltonian system on a Riemannian 3-manifold 
(M,g) satisfying Properties (i) and (ii) in Bertrand's Theorem^ Then the following state- 
ments hold: 

(i) H is of the form ([7]) for some coprime positive integers n, m. 



14 



(ii) The potential V is the intrinsic Kepler or oscillator potential in (M,g). 
(Hi) H is super integr able. More precisely, 

(a) H is geometrically superintegrable in the region ofT*M foliated by bounded orbits. 

(b) There exists an n-fold cover M of M such that the lift of H to M admits a 
generalized Runge-Lenz vector. 

(c) There exists a nontrivial horizontal tensor field in M of rank n which is invariant 
under the flow of H . 

As mentioned in the introduction, this result is of interest both in itself and because of 
the abundant literature devoted to the study of particular cases of this problem in different 
contexts and from various points of view. 
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